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A weak fluctuating magnetic field embedded into a turbulent conducting medium grows expo- 
nentially while its characteristic scale decays. In the interstellar medium and protogalactic plasmas, 
the magnetic Prandtl number is very large, so a broad spectrum of growing magnetic fluctuations is 
excited at small (subviscous) scales. The condition for the onset of nonlinear back reaction depends 
on the structure of the field lines. We study the statistical correlations that are set up in the field 
pattern and show that the magnetic-field lines possess a folding structure, where most of the scale 
decrease is due to the field variation across itself (rapid transverse direction reversals) , while the scale 
of the field variation along itself stays approximately constant. Specifically, we find that, though 
both the magnetic energy and the mean-square curvature of the field lines grow exponentially, the 
field strength and the field- line curvature are anticorrelated, i.e., the curved field is relatively weak, 
while the growing field is relatively flat. The detailed analysis of the statistics of the curvature shows 
that it possesses a stationary limiting distribution with the bulk located at the values of curvature 
comparable to the characteristic wave number of the velocity field and a power tail extending to 
large values of curvature where it is eventually cut off by the resistive regularization. The regions 
of large curvature, therefore, occupy only a small fraction of the total volume of the system. Our 
theoretical results are corroborated by direct numerical simulations. The implication of the folding 
effect is that the advent of the Lorentz back reaction occurs when the magnetic energy approaches 
that of the smallest turbulent eddies. Our results also directly apply to the problem of statistical 
geometry of the material lines in a random flow. 

PACS number(s): 47.27.Gs, 98.35.Eg, 47.65.+a, 05.10.Gg 

I. INTRODUCTION 

It was demonstrated by Batchelor M that a weak magnetic field passively advected by a turbulent velocity field 
would grow, while its characteristic scale would decay. If the magnetic Prandtl number (the ratio of fluid viscosity v and 
magnetic diffusivity r\,Vx = v /rj) is large, there is a broad range of subviscous scales available to magnetic fluctuations, 
but not to fluid motions. This physical situation is realized in such astrophysical environments as the interstellar 
medium and protogalactic plasmas, where Pr ranges between 10 14 and 10 22 , which provides for 7 to 11 decades 
of subviscous range. The weak-field (kinematic) regime is believed to represent the initial stage of the formation 
of the currently observed magnetic fields of galaxies. These fields, which possess a coherent large-scale component 
and whose energies are comparable to the energies of fluid motions of the interstellar medium B-ff , are thought to 
have originated from very weak initial seed fields in the galaxies (or protogalaxies), which have been amplified and 
brought to their current strength and configuration by the dynamo action of the (proto)galactic turbulent plasmas 
(see Refs. |5||| and references therein). Constructing a definitive and quantitative theory of this process remains an 
open problem. This theory must necessarily be a nonlinear one, because the observed fields are not weak. However, 
developing such a nonlinear theory of the magnetic-field evolution will require a thorough understanding of its linear 
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(kinematic) precursor. In fact, this point holds with greater force in view of the recent theoretical and numerical 
advances which suggest that the saturated spectra of magnetic fluctuations are largely determined by the turbulent 
advection processes that drive the kinematic dynamo |J. In this work, we study the geometrical structure of the 
fluctuating small-scale magnetic fields produced by the kinematic stage of the high-Pr dynamo. Our findings will have 
direct bearing on such issues as the condition for the onset of nonlinear effects, the geometry of the field as it enters 
the nonlinear stage of its evolution, and the feasibility of transfering the small-scale magnetic fluctuation energy to 
larger-scale components of the field. 

In an ideal (or highly-conducting) fluid, the magnetic fields are (nearly) frozen into the ambient flow. Therefore, 
besides being important for the astrophysical dynamo as outlined above, studying passive advection of the magnetic 
field is equivalent to studying the statistics of stretching and distortion of material lines by random flows, which is 
a fundamental problem in the theory of turbulence jlO| ]. This subject has attracted considerable attention [10| pl[ . 
Our results on the statistical geometry of magnetic-field lines will have direct applicability in this area. For the sake 
of unity of exposition, we will proceed to develop our theory in the language of the kinematic-dynamo problem and 
relegate the drawing of the parallels with the problem of material-line advection to the end of the discussion section, 
which concludes this work (Sec. 

The mathematical formulation and treatment of the small-scale kinematic-dynamo problem were initiated by 
Kazantsev [^2|. Kulsrud and Anderson (2^] developed a detailed spectral theory of the small-scale magnetic fluc- 
tuations. (A comprehensive exposition of the modern state of the second-order statistical theory of the small-scale 
kinematic dynamo with large Prandtl numbers, as well as the generalization of Kazantsev's and Kulsrud and Ander- 
son's theories to the case of arbitrarily compressible velocity fields, can be found in Ref. |||.) It was established that 
the characteristic scale of the advected magnetic field decreases exponentially fast at a rate comparable to that of 
the field growth. The magnetic spectrum quickly shifts its bulk toward scales extremely small compared to those of 
the velocity field. The decrease of the characteristic scale is checked only by the Ohmic resistive dissipation. Such a 
regime persists as long as the kinematic approximation remains valid. 

It is interesting, and, in fact, necessary for a variety of applications, to inquire what those small-scale fields "look 
like" : do they really tangle into a completely chaotic and fine-scaled web? The most important reason for such 
an inquiry is that it is the structure, not just the strength, of the small-scale magnetic fields that determines the 
conditions for the onset of the nonlinear regime. Indeed, we observe that the Lorentz tension force B • VB only 
involves the parallel gradient of the magnetic field pi] , Heuristically, the nonlinear Lorentz feedback will start to 
play an important role when the Lorentz tension force becomes comparable to inertial terms in the hydrodynamic 
momentum equation, namely, when B 2 ~ (fc„/fc||)pu 2 , where u is the velocity field, k v is the smallest-eddy wave 
number, p is the density of the medium, and k\\ is the characteristic wave number of the magnetic-field variation along 

itself. For chaotically tangled fields, the ratio k v fk\\ can be as small as k v /k n ~ Pr -1 ^ 2 , where k v is the resistive- 
regularization wave number. The kinematic stage of the dynamo will then only produce very weak small-scale fields. 
On the other hand, if fc|| is restricted from growing to be as large as k^, the kinematic dynamo can drive small-scale 
magnetic fluctuations of energies approaching that of the smallest turbulent eddies. Much of the previous work on the 
small-scale dynamo and such issues as ambipolar damping and viscous relaxation of small-scale magnetic fluctuations 
was based on specific assumptions about the magnitude of k 



II ]23| , |25| , |26[ . Understanding the structure of the magnetic 
field, and, in particular, the statistics of the field-line curvature, is also crucial for the study of the effect of the 
Braginskii tensor viscosity p7t] on the small-scale magnetic fields p8|,p9| . 
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FIG. I. Stretching of magnetic-field lines by a linear shear flow. 
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It was suggested by Cowley ||cj on intuitive grounds and later supported by numerical simulations that a 

large-scale advecting field, which is locally just a linear shear flow, could only stretch the magnetic field and make 
it flip direction ever more rapidly in the plane transverse to the field itself [see Fig. |] and Fig. |](a)]. It was argued 
that no appreciable change of the characteristic scales at which the magnetic field varies along itself could, therefore, 
be produced. In other words, the exponential increase of the typical fluctuation wave number k = {k\ + k 2 ) 1 / 2 is 
expected to be due mostly to the increase of k± (rapid transverse direction reversals), while stays approximately 
unchanged, so k± fcii ~ k v . Such folding nature of the small-scale fields is consistent with the predominance 
of volume deformations with greatly disparate spatial dimensions, which is a well known fact in the the theory of 
kinematic dynamo and passive advection pl| , p2[ . It is also, of course, directly related to the extreme flux-cancellation 
property (fine-scale spatial oscillation of field orientation) of the dynamo-generated magnetic fields in maps and chaotic 
flows extensively studied by Ott and coworkers (see review |Q and references therein) and by Cattaneo |34j . 

In this paper, we construct an explicit statistical description of the folding effect in the small-scale kinematic-dynamo 
theory and study the correlations that are set up between the curvature of the magnetic-field lines and the strength 
of the magnetic field. Since we are interested in the geometrical properties of the field, we neglect the resistive effects 
present at extremely small scales and consider the diffusion-free induction equation 

^B = B- Vu-BV-u, (1) 

where 4z = dt + u • V is the full convective derivative, B(i, x) is the passive magnetic field and u(i, x) is the externally 
prescribed velocity field. Let us introduce an auxiliary field F = B • VB, which is, of course, the magnetic-tension 
part of the Lorentz force. It is readily seen that F(i,x) evolves according to the following equation 

— F = F • Vu - 2FV u + BB: VVu -BBVVu. (2) 

dt 

Let us first describe a very simple semiquantitative argument that supports the folding picture. In the incompressible 
case (V • u = 0), we notice that the evolution equation (||) for the Lorentz tension force F is identical to that for 
the magnetic field with the exception of the term BB : VVu, which contains second derivatives of the velocity 
field. Suppose that an initial distribution of the small-scale magnetic fluctuations has been set up in such a way 
that its characteristic parallel and perpendicular wave numbers are comparable and both are much greater than the 
characteristic wave number of the velocity field: k\\ ~ k± 3> fc„. Then the second derivatives of the velocity field can 
be neglected and the mean-square tension force (F 2 ) must grow in the same way as the magnetic energy (B 2 ). For 
the characteristic wave number of the magnetic field we then have 

^ ~ W) K W) < const e ' 12t ' (3) 

where 72 is the growth rate of the magnetic energy (B 2 ) and we have used the obvious fact that (B 4 ) > {B 2 ) 2 . Thus, 
any initial field arrangement where magnetic- field lines are chaotically tangled will decay toward a folding state at the 
rate comparable to the rate of the magnetic energy growth [cf. Fig. |](a)]. 

In order to see how the situation develops when k 2 becomes comparable to k 2 , a more complete analysis of the 
statistics of the magnetic field and the Lorentz tension is required. In Sec. [HJ we we carry out such an analysis exactly 
for the case of incompressible velocity field, and prove that k 2 = (F 2 ) / (B A ) stabilizes at a value ~ k 2 . We then take 



up the question of the evolution of the magnetic curvature, which was recently raised by Malyshkin 35 1. We confirm 
Malyshkin's result on the exponential growth of the mean-square curvature. Most importantly, we find that, while 
the ratio of averages (F 2 ) / (B A ) tends to a constant, the averaged ratio (F 2 j B 4 ) follows the exponential growth of the 
mean-square curvature. This discrepancy implies that the magnetic-field strength and the curvature of the magnetic- 
field lines are very strongly anticorrelated, i.e., the magnetic field is weak wherever the curvature is large, and vice 
versa. The picture of folded magnetic-field lines is manifestly consistent with this property, while that of chaotically 
tangled ones is not. We argue that the large values of curvature in the bends of the folds account for the overall 
growth of the mean-square curvature, even though these bends occupy only a small fraction of the total volume of the 
system. At the end of Sec. [n| we present a simple qualitative description of the folded-magnetic-field-line geometry 
that makes possible the statistical correlations we have found. 



In Sec. Ill, we undertake a more detailed study of the one-point distribution of the magnetic-field-line curvature 
and derive equations for its probability density function (PDF) and all of its moments. This is necessary in order 
to prove the statement of Sec. || that the curvature only grows exponentially in a small fraction of the total volume 
of the system. We discover that, while the moments of the curvature diverge exponentially in time, its distribution 
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tends to a stationary limiting profile whose bulk is concentrated at the values of curvature ~ k v and which has a 
power tail at large values of curvature (the exponent is —13/7 in the 3D incompressible case). We conclude that the 
fraction of the volume where the growth of the curvature takes place tends to zero with time. The limiting values 
of the curvature moments are determined by the resistive regularization at the scales where the magnetic diffusivity 

1/2 

becomes important: k n ~ Pr ' k u . 

Our theoretical results on the folding structure of the magnetic field, the anticorrelation between the field strength 
and the field-line curvature (Sec. |n|), the growth of the curvature moments, and the stationary limiting distribution 
of the curvature (Sec. Ill) are backed up by the numerical evidence based on the 3D incompressible MHD simulations 
by Maron and Cowley ||. The relevant numerical results are reported at the end of each section. The agreement 
between our theory and direct numerical simulations of a realistic MHD environment is quite remarkable, especially 
in view of the idealized character of our modeling assumptions. 

In Sec. IV, we summarize our findings and discuss the implications for the nonlinear dynamo theory. The funda- 
mental consequence of the folding effect (i.e., of the fact that the parallel scale of the small-scale magnetic fields does 
not decay) is that the nonlinear regime sets in only when the magnetic energy becomes comparable to the energy 
of the smallest turbulent eddies. We also explain how our results apply to the problem of statistical geometry of 
material lines in isotropic turbulence and relate our conclusions to the previous work on this subject. Numerical 
results obtained by several authors in this context provide further confirmation of our theory. The fact that the same 
set of basic features of the curvature statistics is found in a number of different approaches and models, many of them 
much more realistic than ours, indicates that these statistics may have a largely universal character. 

The paper also includes three appendices. In Appendix H, we explain the technical details of the derivation of 
the Fokker-Planck equations used in the paper. Appendix [B(is devoted to the study of the structure of the small- 
scale magnetic fields for the case of advecting flows that possess an arbitrary degree of compressibility. We find 
that the folding effect as described above only persists as long as the degree of compressibility of the flow remains 
below a certain critical value. Once this value is exceeded, both the parallel and the perpendicular scales of the 
magnetic-field variation decay exponentially fast (albeit at different rates) into the subviscous scale range and towards 
the resistive scales. If this decay continues until the parallel and the perpendicular scales are equalized, the folding 
pattern is replaced by the tangled one. However, the tangled state is only set up in a small fraction of the total 
volume where the density of the medium is high and where most of the magnetic-field growth takes place. In the 
larger (and less dense) part of the system, the magnetic field stays relatively weak and flat. This new situation 
brought about by compressibility is due to the ability of compressible flows to shrink volumes of the medium with 
frozen-in magnetic-field lines. The structure of the field is determined by the competition between stretching and 
contraction. In Appendix |J the above consideration of the compressibility effects is related to the general theory of 
passive advection in compressible flows developed in Ref. p2[ . 



II. STATISTICS OF LORENTZ TENSION AND MAGNETIC-FIELD-LINE CURVATURE 



In this section, we will restrict our consideration to the case of incompressible velocity field. The evolution equations 
for the magnetic field B(t, x) and the Lorentz tension Y(t, x) in this case are obtained from the equations ([!]) and (^) 
by setting V • u = 0. As is customary in the problems of passive advection |36|,[l2) and kinematic dynamo p^ ], we 
choose the advecting velocity u(t, x) to be a Gaussian white-noise-like random field [the Kazantsev-Kraichnan flow), 
whose statistics are defined by its second-order correlation tensor 

(«*(*, xV'(i',x')) =^-i'K j (x-x'). (4) 

As we will only have to deal with one-point statistical quantities, all the relevant information about the velocity 
correlation properties is contained in the Taylor expansion of k 4 - 7 around the origin j37| 

K iS (V) = *oS ij - i Ka [y 2 S^ + 2ayV ] + ~ n iV 2 [y 2 8^ + 2b y y] + ■ ■ • , (5) 

as y — > 0. In order to ensure incompressibility, we must set a = —\/{d + 1) and b = —2/{d + 3), where d is the 
dimension of space. Our consideration is formally in d dimensions, so that both the two- and the three-dimensional 
cases can be considered in a unified framework. 

The fields B(f , x) and F(t, x) satisfy a a closed system of equations, and, in order to study their statistical properties, 
we derive the Fokker-Planck equation for the joint PDF of B(f, x) and F(t, x) at an arbitrary fixed point x. Due to 
the homogeneity of the problem, this one-point PDF P(t; B, F) is independent of x. A standard derivation procedure 
explained in Appendix leads to the following equation for P: 
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(6) 



The indices following a comma in the subscripts always mean spatial derivatives: : k — d/dx k . K lJ kl and ^\i mn are 

the tensors of second and fourth derivatives, respectively, of the velocity correlator k 4 - 7 (y) taken at y = 0. The 
derivatives with respect to B l and F l in Eq. (^) act rightwards on all terms they multiply. The Einstein convention 
of summing over repeated indices is used throughout. Eq. (||) contains all the one-point statistical information about 
the distribution of B and F and can, therefore, be employed to calculate any individual or mixed averages of these 
quantities. This is done by multiplying Eq. (||) through by the quantity whose average is sought and integrating both 
sides with respect to B and F. The derivatives are removed via integration by parts and an ordinary differential 
equation is established for the desired average, whose time derivative is thereby linked to a linear combination of 
other averages (including itself). The latter averages must in turn be calculated in the same fashion. We will see that 
in many cases of interest, very simple linear equations or closed systems of linear equations emerge. 
Let us start by calculating the mean-square Lorentz tension. We get 

d t (F a )=- lF {F a )+s F {&), (7) 

d t <B 4 H 74 <i? 4 >. (8) 

The expressions for the coefficients 7j?, Sf, 74, as well as for others that will arise in what follows, are collected 
in Table [j]. Note that, in accordance with the simple argument we described in the Introduction, the growth rate jf 
of (F 2 ) is the same as that of the magnetic energy (B 2 ): jp — 72- 
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TABLE I. Coefficients for Sec. |j. The general formulas listed in this table are for the case of arbitrarily compressible flows 
(see Appendix H). The results for the incompressible case considered in Sec. pi are obtained by setting f3 = and £ = 0. 
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Introducing the characteristic parallel wave number of the magnetic fluctuations according to k 2 = (F 2 ) / (B 4 ) , we 
readily find 

Sf \ _( 74 _ TJ? )t , Sf S f k 4 

where k v is the characteristic wave number of the advecting flow. The exponential decay of k 2 was already captured 

in the qualitative argument given in the Introduction [see formula (^)] . The existence of a steady limiting solution is 
due to the presence of the second derivatives of the velocity field in Eq. (|2|). By taking them into account, we have 
thus explicitly proved that k 2 ~k 2 . 

Let us now undertake a slightly more detailed analysis of the magnetic-field structure. The Lorentz tension can be 
decomposed into two orthogonal components 

F = B 2 (b- Vb + b^jj^ =B 2 (K + M), (10) 

where b = B/i3 is the unit vector in the direction of the magnetic field, and V|| = b ■ V. The first term is the magnetic 
curvature vector K = b • Vb, the second term, M = bWnB/B, measures the mirror effect and will, for the sake of 
brevity, be henceforth referred to as the mirror force. Since K _L M, we have (F 2 /B 4 ) = (K 2 ) + (M 2 ). The mean 
squares of both of these quantities can be expressed in terms of mixed averages of F and B: (M 2 ) = ((F • B) 2 /£? 6 ) 
and {K 2 ) = {F 2 /B 4 ) - ((F • B) 2 /B e ), which we proceed to calculate with the aid of Eq. (§) 

d t (K 2 ) = 1K (K 2 ) + S K , (11) 
d t (M 2 } = - 1M {M 2 ) + 1MK {K 2 ) + S M (12) 

(see Table | for the values of the coefficients). The exact solution of Eq. ( O) is 

- ( iTs2\ ln \ , Sk \ „iKt Sk 



(K 2 )(t) = <A 2 }(0) + ^ e w* - ^, (13) 
V Ik) Ik 

so the magnetic curvature grows exponentially (even if it is initially zero). It is instructive to express its growth rate 
in terms of the growth rate of the magnetic energy (B 2 ). In three dimensions, this gives — (16/5)72/2, which 
agrees with the result Malyshkin |55j obtained by a direct calculation of (A" 2 } in the spirit of the Kulsrud- Anderson 
theory [ p3| . We also see that the mean-square mirror force [Eq. ([l2|)1 is not an independently interesting quantity: 
after a transient initial time, it is reduced to "mirror" the evolution of the mean-square curvature: 

(M 2 )(t)~^^(K 2 )(t), t-00. (14) 

IK + 1M 

Thus, we have established that, while the ratio of the averages (F 2 ) / (B A ) tends to a constant value ~ /c 2 , the 
averaged ratio (F 2 / B 4 ) ~ (A 2 ) ~ e lKt grows exponentially. Since both of these quantities have the dimension and 
the intuitive meaning of some characteristic parallel wave numbers, the question inevitably arises as to the physical 
interpretation of such drastic dependence on the relative order of the averaging and the normalization with respect to 
the magnetic-field strength. This dependence clearly indicates that there exists a very strong anticorrelation between 
the strength of the magnetic field and the curvature of the magnetic-field lines. Namely, while both the mean-square 
curvature and all moments of B grow exponentially, the magnetic fields are configured in such a way that the magnetic 
field is very weak wherever its curvature is large, and vice versa. No such arrangement would be possible if the field 
were chaotically tangled everywhere. Indeed, a tangled state of this sort would imply that the absolute values of the 
curvature were everywhere comparably large and growing. But then, in order to compensate for the growth of the 
curvature, the growth of B 4 would have to be partially or fully suppressed compared to that mandated by Eq. (JsJ) . 

On the other hand, it is easily envisioned how the strong anticorrelation between B and A can be realized in the 
folding picture. While the curvature is quite small and magnetic field grows in most of the volume, which is occupied 
by the folds, the situation is reversed in the small part of the volume where magnetic-field lines bend and reverse 
direction: the curvature there is very large and magnetic field weak. Fig. illustrates the typical geometry of the 
folding magnetic- field lines. This picture is in 2D, but can also be interpreted as a cross section of a 3D geometry of a 
sheetlike configuration. Flux conservation (J B • dS = 0) implies -Bbcnd/Sfoid ~ £±/h, where £± is the characteristic 
scale of magnetic-field variation across itself in the folding region and £b is the characteristic size of the bend. The 



G 



velocity shear that produces (or "sharpens" ) the bend acts in such a way that £± is decreased while 4 is amplified, 
so 4 > £±, whence B hcnd < -Bfoid- 

B=0 




FIG. 2. The geometry of the folding field lines in the vicinity of the bend. The dashed lines correspond to the surfaces on 
which the magnetic field vanishes. The shaded area is the cross section of the volume that can be used for the flux-conservation 
estimate -E>b cn d/-Bfoid ~ i±/^b- All the flux is through the surfaces whose cross sections are depicted by the bold lines. 




(a) Magnetic field B (b) Curvature K = |b • Vb| 



FIG. 3. Instanteneous magnetic-field configuration in the kinematic regime (numerical results). These are greyscale plots of 
2D cross sections of the 3D snapshots of (a) the magnetic-field strength and (b) the absolute value of the field-line curvature. 
These plots are from the same simulation as Fig. ^(b). The field-strength and the curvature snapshots [plots (a) and (b)] arc 
taken at the same moment t — 8.2 and at the same cross section. Darker regions correspond to larger values of the fields. In the 
plot (a), (B) ~ 0.003, (B 2 ) 1/2 ~ 0.004, (B i ) 1/i ~ 0.006, and the maximum value of B throughout the system is ~ 0.025. The 
regions that are pitch black in the plot encompass fields stronger than 0.01. All of these values correspond to magnetic-field 
energies well below the nonlinear-saturation threshold. The specific units of the field strength are, of course, of no consequence 
here. In the plot (b), (K) ~ 50, {K 2 ) 1/2 ~ 70, {K 4 ) 1/4 ~ 110, and the maximum value of K is ~ 520. The pitch black regions 
of the plot correspond to curvatures larger than 400. The curvature has units of inverse length, based on the box size 1. 
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It must be recognized, however, that the presence of an anticorrelation between the magnetic-field strength and 
the magnetic-field-line curvature does not in itself prove that the volume where the growth of the curvature occurs 
constitutes only a small fraction of the total volume of the system. Indeed, examples of magnetic fields can be 
constructed that possess such an anticorrelation and where, at the same time, the mean-square curvature grows in 
any arbitrary fraction of the total volume that can be specified beforeh and. Further study of the curvature statistics 
is, therefore, required to settle this issue. This will be carried out in Sec. Ill, where the smallness of the volume where 
the curvature grows is confirmed. 

Finally, let us reiterate that the presence of the folding structure has found repeated confirmation by numerical 
evidence. Most recently, folding was extensively studied in 2D and 3D numerical simulations of the small-scale dynamo 
effect in a viscosity-dominated MHD model of Kinney et al. , and in 3D forced-MHD simulations of Maron and 
Cowley H . Here we present the numerical results that are based on the latter work and directly relate to the theory 
developed in this section. All numerical results presented in this paper derive from a 128 3 spectral forced-MHD 
code written by J. Maron and described in detail in Refs. |3^|J. The external forcing is on the system-size scale 
and S correlated in time. In the simulations quoted in this paper, the hydrodynamic Reynolds numbers are quite 
small, so the advecting fluid flows are essentially determined by the balance of the forcing and the viscous dissipation. 
However, this is not really a handicap, as the purpose of the numerical results presented here is to illustrate the 
kinematic-dynamo properties at subviscous scales. More discussion of this issue can be found in Refs. RS. 
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(a) Decay of (F 2 )/(B 4 ) (log 10 plot) 



(b) (K 2 ) and (F 2 )/(B 4 ) (log 10 plot) 



FIG. 4. Anticorrelation between the magnetic-field strength and the field-line curvature and growth of the mean-square 
curvature (numerical results), (a) Time evolution of (F 2 )/(B 4 ) in two simulations where the initial magnetic field is uniformly 
tangled at subviscous scales. The hollow dots correspond to the simulation with Pr = 100, k v ~ 25; the filled dots correspond 
to the simulation with Pr = 2500, k v ~ 5. (b) Time evolution of (F 2 )/(B 4 ) (hollow dots) and {K 2 ) (filled dots) in a simulation 
with the initially flat magnetic field varying transversely at the velocity scales (fcii = 0, k± ~ k v ). In this simulation, Pr = 2500, 
kv ^ 5, k-q ~ 250. The ratio (F 2 )/(B 4 ) again stabilizes at a value ~ a few times k 2 . In both plots, the quantities plotted 
have units of inverse length. These are based on the box size 1. Time is measured in units of the smallest-eddy turnover 
time (|V x ul 2 )" 1 ^ 2 (in units based on box size 1 and forcing power 1, this quantity is ~ 0.22). 

In Fig. ||, we give the greyscale plots of the magnetic-field strength and the absolute value of the magnetic curva- 
ture corresponding to a typical instantaneous magnetic-field configuration observed in a simulated MHD evironmcnt 
during the kinematic stage of the small-scale dynamo. The folding pattern strikingly similar to the one described 
above is clearly in evidence (cf. Fig. [2]). The anticorrelation between the field strength and the field- line curvature, 
as well as the intermittent nature of the distribution of both (cf. Sec. ill) are also manifest. Fig. ||(a) shows how the 
ratio k 2 = (F 2 ) j (B A ) adjusts to a stationary value ~ k 2 from an initial state where the field is chaotically tangled 
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at subviscous scales. We show results of two simulations with such an (artificial) initial field and different values 
of k u . In both cases, exponential decay of fc|(t) toward stationary values ~ a few times k 2 is observed, which corrob- 
orates Eq. (||). Fig. [|(b) portrays the time evolution of the ratio (F 2 ) / ' (B A ) and of the mean-square curvature (K 2 ) 
in a simulation that starts with the magnetic field concentrated at the velocity scales. The ratio (F 2 ) / (B A ) again 
stabilizes at a value ~ a few times k 2 as predicted by our solution (j^). The exponential growth of the mean-square 
curvature (K 2 ) proceeds in accordance with our solution ([l3]) until it is checked by the resistive regularization at a 
stationary value (K 2 ) ~ k 2 . While our theory was constructed for the diffusion-free regime and, therefore, did not 
include this effect, the resistive saturation of the curvature is naturally an expected outcome (see Sec. [II for more 
discussion of this issue). 



III. DISTRIBUTION OF MAGNETIC-FIELD-LINE CURVATURE AND MAGNETIC-FIELD STRENGTH 

In the previous section, we indicated the need for a study of the curvature statistics that would go beyond the 
evolution of the mean square. In this section, we fulfill this program and delve deeper into the detailed properties of 
the distribution of the magnetic field and its curvature. 

The Fokker-Planck equation for the one-point PDF of the magnetic-field-line curvature K = b • Vb is most conve- 
niently derived on the basis of the folowing coupled evolution equations for K and the magnetic-field direction b: 

K = K • (Vu) • (i - bb) - bKb : Vu - 2Kbb : Vu + bb : (VVu) • (I - bb), (15) 

^b = b-(vu)-(i-bb), (i6) 

where I is the unit dyadic and colons denote double dot products executed according to Kb : Vu = K • (b • Vu), etc. 
Both of the above equations are direct corollaries of the induction equation (0) . It is easy to see that these equations 
respect the conservation laws b — 1 and b • K = 0. Note that, in this section, we work with arbitrarily compressible 
velocity fields, so u is not required to be divergence-free. It will be seen, however, that none of the essential features 
of the curvature statistics are affected by the compressibility. 

The averaging procedure that leads to the Fokker-Planck equation for the joint PDF P(t; K, b) does not involve 
any nonstandard steps and is fully analogous to that used to derive the Fokker-Planck equation (||) (see Appendix |a|). 
The result is 

d t P = -I K ij kl (-Si + -^b k + -^—K k - -^-b r b k b l - - -^—WK k V - -^—Wb k K' 

2 > kl \ 1 db l dK l db r dK r dK r dK r 

dbi dKi db s dK s dK s dK s 

1 K % ( JL_b k b m - -^—b r b k b m b^\ ( -^—b l b n - -^—b s b l b n b>\ P. ( 1.7) 

2 ' fcimn V dK l dK' r / V dK J dK s 1 



A major simplification of this equation becomes possible if one recalls that the joint distribution P(t; K, b) is subject 
to two constraints: |b| = 1 and b • K = 0. Also taking into account the spatial isotropy of the problem, we conclude 
that the following factorization holds: 

P(t; K, b) = ^(|b| 2 - l)<5(b • K)P K {t; K). (18) 
The function Px{t] K) is then found to satisfy the following reduced Fokker-Planck equation: 



dtPK = 2(5Ti) K2 



5d - i + k 2 P'^ + (lid 2 -6d+l+ 2{7d d 2) (3 ) K1>' K 



where primes denote partial derivatives with respect to K , the compressibility parameter (3 = d[l + [d + l)a] is 
nonnegative and vanishes in the incompressible case, and «2, K4, and a are coefficients of the small-scale expansion (^) 
of the velocity correlator. Note that the distribution of the curvature is independent of the second compressibility 







parameter b. The normalization rule for Pff(i; K) follows from the normalization of the original PDF Pit; K, b) and 
from the factorization ©: (l/2)S d S d -i J °° dK K d - 2 P K (t; K) = 1, where S d = 2n d / 2 /T(d/2) is the area of a unit 
sphere in d dimensions. Absorbing the geometrical prefactor into PK(t;K), we conclude that the true PDF (in the 
sense that it induces a measure on the volume of the system and integrates to unity) is K d ~ 2 PK(t; K). Note that, 
since the curvature vector must always remain perpendicular to the direction of the magnetic field, the curvature 
distribution is effectively restricted to d — 1 dimensions. 

It is now straightforward to establish the set of evolution equations for the even moments of the curvature 



d t (K 2n ) = 



2(5d-l) 



1 



-d 



2 

d V d 



Gn 



1 



-1/3 



riK 2 (K 2n ) + 6{d + 2n- 3)nK i (K 2(n - 1) ) , n > 1 



(20) 



For n — 1, Eq. (|2(]) reproduces the results for the mean-square curvature that were obtained in Sec. [n| and Appendix [b| 
[see Eq. (^) and Table The higher moments of the curvature are coupled to the lower ones in a recursive fashion, 
but also have their own growth rates that increase quadratically with n. This latter kind of intcrmittency is very 
similar to that encountered in earlier studies of the statistics of the magnetic- field strength J3l],|32j . For the sake of 
comparison, let us list here the Fokker-Planck equation that determines the PDF B d ~ 1 PB(t; B) of the magnetic-field 
strength B and the evolution equation for its moments (B 2n ) = Sd / °° dB B d ~ 1+2n Ps{t; B), 



d t P B 



1 d 

— Ko 

2 



1 



d+1 

d t (B 2n ) = 



(1 + p) B 2 P'^ + (d+ 1)(1 + 2(5)BP' B + d{d + 1)/3P B 
d-1 



- [2n + d+{2n-l)f3]nK 2 (B 2n ). 



(21) 
(22) 



The primes in Eq. (|T]) denote derivatives with respect to B. We note that Eq. (||) is a particular case of Eq. ^2%). 
Direct derivation of the above equations by averaging the induction equation ([j]) is quite standard. Details can be 
found in Ref. J|s|. Eq. (|2l]) can also be obtained by integrating out the F l dependence in Eq. (B2) and using the 
spatial isotropy of the magnetic-field distribution. Eq. ( p2] ) is a direct consequence of Eq. (pi]). 

We now turn to the main objective of this section, namely, estimating the fraction of the total volume of the system 
where the curvature growth occurs. In Eq. (pj|), denote by D, E, and T the coefficients in front of K 2 P^, KP' K , and Pk, 
respectively. Now rescale time and curvature according to Dt => t and K/K* K, where = (3^4/ D) 1 ! 2 ~ k v 
(recall that fc„ is the characteristic wave number of the advecting velocity field). We can now rewrite Eq. (19) in the 
following nondimensionalized form 



d t P K - (1 



K 2 )P'^ 



aK 



d-2 
K 



P' 



{d-l){<J~d)P K , 



where we have used the fact that T = (d— 1)(S — Dd) and denoted 

S lid 2 -M+l + 2{7d-2)[3/d 



D 



5d-l + 6(3/d 



(23) 



(24) 



Besides the dimension of space, the only essential parameter of the curvature distribution is a, which changes with d 
and the degree of compressibility. The correct boundary conditions for Eq. ( ^3| ) follow from the normalizability 
requirement dK K d - 2 P K (t; K) < 00, 



[K d - 2 P' K (t; K)] K=0 = 0, [K d P' K (t; K) + (a - ^K^Pk^ K)] = 



(25) 



Let us study the evolution of the curvature statistics from an initial setting where the curvature is zero every- 
where: K d ~ 2 PK{t — 0, K ) oc 5(K). While such a 8-\ike initial distribution is, of course, highly artificial, mathemati- 
cally it is not an anomalous case since, as we have seen [Eq. (p0[)], the moments of the curvature would grow even from 
such an initial state. Two distinct asymptotic regimes can be identified in the evolution of the curvature distribution. 

Small-curvature regime. For small values of curvature K <^L 1 (i.e., for the dimensional curvature much smaller 
then K* ~ k v ), Eq. (|23|) reduces to what mathematically is a heat equation in d — 1 dimensions with radial symmetry 



d t P K = P'k 



d-2 
K 



P'k- 



(26) 



The solution is a heat profile spreading out from the origin 
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e -K 2 /4t 

P K (t;K) = const t(d _ 1)/2 - (27) 



Multiplying the solution ( |27j ) by K d ~ 2 , we find the peak of the PDF at K pca ^ = \/2{d — 2)t, i.e., it remains at K = 
in 2D and shifts towards larger K in 3D. In either case, the excitation eventually spreads over towards larger K ~ 1, 
where the small-curvature asymptotic regime breaks down. 

Large-curvature regime. At large values of curvature K ^> 1, the asymptotic form of Eq. (f23|) is 



= K 2 P£ + aKP' K + (d - 1)(<t - d)P*. (28) 

In logarithmic variables, this is a ID diffusion equation with the drift velocity a — 1 and with an overall growth 
rate (d — 1)(<t — d). The corresponding Green's function is 

e ( d -i)( g - d )( t - t0) / [ m (jf/if ) + ( ff -l)(t-t )] 2 ^ 

G K {t-t ;K,K ) = cxp -i — — — I ■ 29 

K 0y /4Tr(t-t ) \ 4(*-t ) y 

Thus, the curvature distribution develops a lognormal tail. This clearly accounts for the intermittency we have de- 
tected in the evolution of the curvature moments [Eq. (|20|)]. Multiplying the Green's function (p9|) by K d ~ 2 , it is not 
hard to see that the peak of the excitation propagates according to K pca k = K$ exp [(2d — 3 — <r)(t — to)] . Substituting 
the value of a [formula (pi[)], we see that 2d — 3 — a < in both two and three dimensions, so the peak, in fact, 
propagates backwards towards smaller values of curvature. 

The conclusion from this simple asymptotic analysis is that, after an initial transient time, the curvature PDF 
should assume the form where its bulk is concentrated at the values of (dimensional) curvature smaller or comparable 
to ~ k v and a lognormally decaying tail is formed at K 3> K*. The global maximum of the PDF is located 
at K — in 2D and at some K ~ in 3D. Thus, in most of the volume of the system, the values of the curvature 
should not greatly exceed K*. The growth of the moments of the curvature is, on the other hand, mostly due to 
the lognormal tail of the distribution. Indeed, multiplying the solutions (^7]) and ( p9| ) by J{ d '~ 2+2n , we see that the 
relative importance of the small-curvature region decreases, while that of the lognormal tail increases. The peak of 
the function K d ~ 2+2n Pk{K) for n > 1 always propagates in the forward direction. 

Of course, once the solution of Eq. ( |2^ ) has cleared the region of validity of the small-curvature asymptotic regime, a 
complicated process of probability redistribution is set up. As time passes, the lognormal tail gains more weight, while 
the heat profile at small K spreads out. The precise nature of the evolution of the PDF is decided by the interaction 
between the small-curvature (radial-heat) and large-curvature (lognormal) regimes in the crossover region K ~ K*. 
This interaction can affect the entire PDF. We can gain more insight into what happens by observing that Eq. (]23|) has 
a stationary solution. Indeed, let us write Eq. ( p3| ) in the following explicitly conservative form: 

Q t P K = J^2^ Rd ~ 2 [(! + K ") P K + (° ~ d ) KP K\ ■ ( 3 °) 



Setting the left-hand side to zero, integrating twice, and making use of the boundary conditions (25) to eliminate one 
of the constants of integration, we find the following stationary limiting PDF: 

K d - 2 P^\K) = const (1 + ^L)/2 - ( 31 ) 



This PDF satisfies the boundary conditions (p5[), is properly normalizable, and has a power tail ~ /^ _ [ cr_2 ( d_1 )]. 
The values of the exponent a — 2(d — 1) for the incompressible and irrotational cases in two and three dimensions, 
along with the values of other relevant parameters, are collected in Table |l[ The curvature distribution can be seen 
to converge (in the mean-square sense) to the stationary profile (^TJ) if we represent the time-dependent solutions 
of Eq. ( |30| ) in the form Pi((t;K) — C(t; K)P^ t \K) and notice that the prefactor C(t;K) tends to a constant, viz., 
(9 f (C 2 ) = — 2((1 + K 2 )(dC/dK) 2 ^, where the averages are with respect to the stationary distribution (|3l|). 

Fig. H shows the results of numerical solution of Eq. (|3(]) in two and three dimensions for the case of incompressible 
velocity field. The results for the case of irrotational velocity field are very similar in form (with slightly different 
power laws at large K: see Table |ll|). Collapse of the curvature PDF onto the stationary profile ([H]) is very quick 
and proceeds in essentially the same fashion in both incompressible and irrotational cases (see discussion at the end 
of this section). 
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InK 




(c) d = 3, KP K (t; K) (d) d = 3, KPr:^; K), ln/ln plot 

FIG. 5. The results of a numerical solution of Eq. in two and three dimensions for the case of incompressible ve- 

locity field. The numerical solution was initialized with the Gaussian heat profile ( |27| ) corresponding to t = 0.2. Time is 
measured in the units of D^ 1 , curvature in the units of K*. We plot the evolving PDF K d ~ 2 PK(t;K) (normalized to 1) at 
times t = 0.2,0.5,1,2,3,4,5,6,7,8,9. The PDF that corresponds to the earliest time is the one with the highest peak and 
the steepest decay at large K. At later times, the peak of the PDF descends, while the tail becomes thicker (lognormal with 
increasing variance and eventually powerlike). The log/log plots (b) and (d) illustrate how the power tail is formed. For the 
sake of reference, we have also plotted the slopes correspoding to K~^ a ~ 2( - d ~ 1 ^ (see Table Inj). 



Velocity Field 


Dimension 


a 




D/K3 


E/K2 


F/K2 


a 


a -2(d-l) 


Incompressible 


d = 3 


-1/4 





7/4 


41/4 


10 


41/7 


13/7 




d= 2 


-1/3 





3/2 


11/2 


5/2 


11/3 


5/3 


Irrotational 


d = 3 


1 


15 


11/2 


34 


35 


68/11 


24/11 




d= 2 


1 


8 


11/2 


43/2 


21/2 


43/11 


21/11 



TABLE II. The coefficients of Eq. @ and Eq. @. 



12 



50 

K 



(a) d = 3, 4TT 2 KP K (t; K) (b) d = 3, 4^ 2 A^P K (t; X), log 10 /logi plot 

FIG. 6. Curvature PDF from 3D incompressible MHD simulations. This is the same simulation as was used 
in Fig. §(b): Pr = 2500, k v ~ 5, k v ~ 250, flat initial field. We plot the PDF An 2 KP K (t; K) (normalized to 1) at 
times f = 0.9, 1.8, 2.7, 3.6, 4.5, 5.5, 6.4, 7.3. The curvature has units of inverse length, based on box size 1. The time is measured 
in the units of (|V x u| 2 ) _1//2 (the smallest-eddy turnover time). The PDF that corresponds to the earliest time is the one 
with the highest peak and the steepest decay at larger K. At later times, the peak becomes lower, the decay gentler, (a) is the 
linear plot, (b) is the log/log plot. The dashed line represents the slope corresponding to A~ 13//7 ' . 

An important feature of the stationary PDF (|3l|) is that all the moments (K 2n ) diverge. In the language of 
physical reality, this means that the limiting values of the curvature moments are essentially determined by the 



resistive regularization, which must cut off the power tail of the PDF ( |31[) at the scale where magnetic diffusivity 

1 /2 

becomes important: k n ~ Pr ' k v . This is, of course, hardly surprising because curvature is just a measure of the 
inverse scale of the magnetic fluctuations and cannot exceed the resistive scale. In view of these findings, the growth 
rates for the curvature moments that have been obtained in this and the preceding sections, should be interpreted as 
describing the evolution of the moments while the lognormal tail of the evolving distribution spreads and thickens. 
The power tail of the stationary limiting distribution |3l]) forms the envelope inside which this process takes place. 
In the diffusion-free regime, the stationary distribution itself is attained at t — > oo with the moments diverging 
exponentially in time. We should like to observe here that a PDF such as we have obtained, with a power tail and 
divergent moments, is indicative of a fractal nature of the distribution. It must be clear that the presence and the 
particular form of the small-scale regularization may affect the global shape of the curvature distribution. Since we 
work in the diffusion- free limit, our theoretical results only apply to the period in the evolution of the magnetic 
fluctuations before the small-scale cut off is reached. In astrophysical plasmas with very large Pr, this corresponds to 
an appreciable length of time. In fact, current estimates pertaining to the (proto) galactic dynamo suggest that the 
kinematic approximation may well break down before the resistive scales become important H . 

As regards the distribution of the curvature over the volume of the system, the fraction of the volume where the 
curvature exceeds any given value Kq is easily seen to be V(K > Kq) — J K dK K d ~ 2 PK{K). The existence of a 
stationary distribution implies that this quantity tends to a constant that depends on the value chosen. Since the 
bulk of the distribution remains at the values of the (dimensional) curvature comparable to AT* ~ fc„, the value 
of V(K > Kq) for Kq 3> AT* will be small. For example, we can use the stationary distribution (|3|) to estimate that, 
in the 3D incompressible case, the fraction of the volume where the curvature is more than ten times larger than 
does not exceed 14%, while the fraction of the volume where the curvature is larger than 100 A'* is no more than 2%. 
It is, of course, quite clear that the fraction of the volume where exponential (or any other kind of) growth of the 
curvature occurs tends to zero with time. 
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Now let us compare the properties of the curvature distribution we have just described with the properties of 
the PDF of the magnetic-field strength determined by Eq. (|2l|). Clearly, the Green's function Gsit — to; B,Bo) for 
this equation is everywhere lognormal and analogous in form to the function Gx{t — t ; K, Kq). With time rescaled 
according to tK 2 (d — l)/2(d + 1) => t, we have 

w-usb,*) :" d+im -' ^HT?-*"! (32) 

BoVS(l +P)(t-to) \ 4(l+/))((-(„) J 

Multiplying Gs(t — to] B, Bo) by £? d_1 , we find that the peak of the excitation is at £?p C ak = Bq cxp[(d — 2 — 3/3) (t — to)] ■ 
In the incompressible case, the peak is stationary in 2D and propagates forward in 3D. For all compressible flows 
in 2D and those with a > —2/9 in 3D, the direction of propagation is reversed. On the other hand, the peak of the 
function B d ~ 1+2n Pb(UB) for n > 1 always propagates in the forward direction, which accounts for the growth of the 
moments of B [Eq. (|22|)1 . 

Thus, we see that the weakening of the anticorrelation between field strength and curvature in the compressible 
flows (Appendix |b]) is due not to any essential change in the properties of the curvature distribution, which is quite 
insensitive to the variation of the degree of compressibility of the flow, but rather to the fact that the magnetic-field 
strength itself now tends to only grow in a decreasing fraction of the total volume of the system. In such a case, both 
the curvature and the magnetic field remain relatively weak in most of the volume. An essential difference in their 
statistics is that, unlike the curvature, the magnetic field does not possess a stationary limiting distribution. Indeed, 
due to the scale-invariant nature of the Fokker-Planck equation (|2l]), such a distribution would have to be a global 
power law and hence could not be normalizable. 

In conclusion, we check the main results obtained in this section against the numerical evidence supplied by the 
3D incompressible MHD simulations of Maron and Cowley In Fig. [| we present the evolution of the curvature 
distribution observed through a sequence of times that correspond to the kinematic and diffusion-free stage of these 
simulations (at later times, the PDF is affected by the resistive regularization and then by the nonlinear effects). 
The collapse of the curvature PDF onto a stationary profile is manifest. The log/log plot [Fig. ||(b)] confirms the 
emergence of the power tail ~ X~ 13 / 7 . These results agree very well with our theoretical predictions presented in this 
section. 



IV. SUMMARY AND DISCUSSION 



Let us now summarize the main physical points we have pursued in this work and discuss the implications for the 
nonlinear dynamo theory. In this section, we will only discuss the case of incompressible flows as the most relevant 
in the astrophysical context we have in mind. The effects of compressibility have been given ample attention on both 
quantitative and qualitative level in Appendices |^ and |^ and in Sec. III. 

In the astrophysical environments that have large magnetic Prandtl numbers and, therefore, possess a wide range of 
(subviscous) scales available to the magnetic, but not hydrodynamic, fluctuations, the small-scale kinematic dynamo 
is driven by the velocity field that locally looks like a linear shear. The volume deformations produced by this field 
lead to exponentially fast stretching and folding of the magnetic-field lines into a structure characterized by very rapid 
transverse variation of the field, which flips its direction at scales ultimately bounded from below only by the resistive 
length. However, the field lines remain largely unbent up to the scale of the advecting flow. 

Both numerically and analytically, we have established that the curvature of the magnetic-field lines and the field 
strength are anticorrelated, i.e., the growth of the field (dynamo) mostly occurs in the regions of flat field while the 
sharply bent fields remain relatively weak. This situation is quickly restored even if the field is artificially scrambled 
into a chaotically tangled state. Moreover, in the three-dimensional incompressible flows, it is the flat growing fields 
that occupy most of the volume of the system. Accordingly, the field-line curvature remains comparable to the inverse 
velocity scale in most of the volume, though its distribution is intermittent and all of its moments grow exponentially 
on account of the small regions of strongly bent (but weak) fields. 

In the diffusion-free approximation, i.e., in the regime where the magnetic excitaion has not reached the resistive 
scale, the growth of the curvature moments is unbounded and the curvature distribution tends to a stationary limiting 
power-like profile with divergent moments. If the resistive cut off is felt while the magnetic field is still weak enough to 
satisfy the kinematic assumption, the curvature moments saturate at the resistive scale and the precise global shape 
of the curvature distribution may be modified. However, the main features of the folding structure described above 
(the anticorrelation between the curvature and the field strength, the smallness of the volume where the field is bent) 
survive because they result from the large-scale geometric properties of the advection rather than from the particular 
form of the small-scale regularization. 
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Let us remark that the reason for some of the statistical quantities considered in this paper achieving steady- 
state values even within the confines of the diffusion- free kinematic approximation ((F 2 ) / (B 4 ) — > const, stationary 
curvature PDF) is that the second derivatives of the advecting velocity field appear in the corresponding dynamic 
evolution equations [see Eq. (||) and Eq. (|l||)]. While passive fields such as B only feel the linear component of 
the ambient velocity field and, therefore, have scale-independent nonstationary distributions, the statistics of their 
gradients involve an additional scale-dependent parameter K4/K2 ~ k 2 . In other words, whereas the magnetic field 
only knows that it is advected by a large-scale flow, the statistics of the magnetic-field gradients specifically depend 
on the actual scale size of this flow and would not be fully captured in a theory where the velocity field were assumed 
to be linear. Physically, the first derivatives of the velocity cause the stretching of the field lines (random shear), 
while the second derivatives are responsible for the bending. 

The theoretical results presented in this paper have been derived for the Kazantsev-Kraichnan model velocity 
field (|J), which is Gaussian and S correlated in time. In view of the highly artificial nature of this field, the question 
naturally arises as to the validity of such results in application to physical, or even to numerical, realities. Indeed, the 
velocity field that arises in the Kolmogorov turbulence is neither Gaussian nor S correlated in time: its intermittent 
character is well known and its correlation time is naturally estimated to be of the same order as its eddy-turnover 
time. However, the Kazantsev-Kraichnan model, while by no means a controlled approximation ]39|] , appears to 
correctly capture most of the physics of the passive advection on at least a semiquantitative level and has survived 
a number of reality tests, numerical and, in the case of scalar turbulence, also experimental (see, e.g., the 

review f40(| and references therein). Indeed, as was reported in the preceding sections, our results are in a very good 
agreement with numerical simulations, where the velocity field derives from the forced Navier-Stokes equation and 
has a realistic correlation time. 

The primary motivation of this study of the structure of the magnetic field was its crucial importance for the 
understanding of how the nonlinear effects set in. We reserve the detailed qualitative and quantitative discussion of 
this issue for an upcoming publication . Here we restrict ourselves to mentioning the most immediate consequence 
that the folding nature of the small-scale field in the kinematic regime has for the onset of the nonlinearity. The 
Lorentz-feedback term in the MHD momentum equation is proportional to the Lorentz tension force F = B • VB [B4| . 
This quantity is quadratic in the magnetic-field strength and involves the parallel gradient of the field. The overall 
effect of the correlations that produce the folding structure is to fix the effective value of this parallel gradient at 
approximately the inverse velocity scale fc„. Thus, the condition for the nonlinearity to become important is the 
growth of the magnetic energy to values comparable to the energy of the smallest turbulent eddies, rather than to much 
smaller values at which the Lorentz tension of a chaotically tangled field would start balancing the inertial terms 
in the momentum equation. Any prospects for producing magnetic fluctuations at larger scales depend on whether 
there exists a nonlinear mechanism for unwinding the folded structure that the nonlinear regime inherits from the 
kinematic one. For further discussion of this subject the reader is referred to Refs. 

Finally, as was promised in the Introduction, let us discuss the relation of our results to the fundamental turbulence 
problem of material-line advection. In an ideally conducting fluid, the behavior of the magnetic-field lines and that of 
the material lines are, of course, identical, so the diffusion-free kinematic-dynamo problem can be recast as a problem 
of stretching of the material lines by the ambient flow. The pioneering work on this subject is due to Batchelor Jl0| , 
who realized that, on the average, turbulent motions lead to exponential elongation of material line and surface 
elements. His results were extended and, in part, made rigorous by a number of authors [Tl|-[l5|]. Starting with the 
work of Pope |l(f | , much attention was focused on the statistics of the extrinsic principal curvature of material surface 
elements in turbulent flows JT~6|— |l9[| and of the curvature of material lines in both turbulent Jl8|,|l4|,[l5| and deterministic 



but chaotic |20| flows. Our results on the statistics of magnetic-field- line curvature are subject to direct comparison 
with these earlier studies. Indeed, our Eq. ( |Hf ) for the vector curvature K = b • Vb of the magnetic field line leads 
to the following equation for K = |K|: 

— K = -(2bb : Vu - fin : Vu)K + bb : (VVu) • n, (33) 
at 

where n = JZ/K is the unit normal to the field line. Eq. ([33]) is the same as the equation derived by Drummond and 
Munch p8| for the curvature of the material line elements (in the usual geometric definition) and is very similar (though 
not identical) to Pope's |l6| equation for the extrinsic principal curvature of the material surface elements. Formal 
similarity between the curvature equations for material lines and surfaces led Drummond and Munch []l8| to conjecture 
that, in isotropic turbulence, the general features of the curvature statistics for these objects would also be alike [ pT| . 
All of the extant numerical evidence supports this conjecture at least on the qualitative level [ p7Ul^JT^ , [l5| , ^o|] . The 
curvature statistics also appear to be largely insensitive to the type of flows considered. Thus, numerical simulations 
involving 3D forced jl7| and 2D decaying pil N avier-Stokes turbulence, Kraichnan's random flow model |p^| , |l5| , 
2D and 3D deterministic but chaotic flows pU|, and, finally, our own 3D forced-MHD simulations and theoretical 
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results based on the Kazantsev-Kraichan velocity field (both compressible and incompressible), consistently reveal 
the same set of properties of the curvature distribution. The unbounded exponential growth of the mean-square 
curvature was first observed in the numerical studies of Pope and coworkers flTj], who also found that, unlike the 
moments of the curvature itself, the moments of its logarithm tended to time- independent asymptotic values. The 
emergence of a stationary curvature PDF with a power tail was reported. Drummond [jl5| confirmed these results in 
his numerical model and also offered a qualitative theoretical argument that related the existence of the power tail 
of the curvature PDF to the competition between the stretching action of the first spatial derivatives of the velocity 
field and the bending effect of its second derivatives. Ishihara and Kaneda (IJ] attempted to determine the power law 
(in 2D) by looking for a critical index p c such that the curvature moments of orders higher than p c would diverge. 
For a model velocity field essentially equivalent to the Kazantsev-Kraichnan flow, they deduced p c = 2/3, which 
was consistent with their numerical results for a realistic (decaying) 2D turbulence. All of the above is in perfect 
agreement with the results of Sec. [II (in particular, Ishihara and Kaneda's |lj] critical index exactly corresponds 
to our —5/3 power tail in d = 2). Drummond and Munch ]l8[ ] also (numerically) measured the correlation between 
curvature and stretching and found it negative, as did we in our theory and simulations (Sec. |J). Such anticorrelation 
was also noticed by Boozer and coworkers |2l|, who, in their theory of finite-time Lyapunov exponents for chaotic 
flows, found that the local Lyapunov exponent of the flow was strongly suppressed in the regions of high curvature. 
Finally, results very similar to those surveyed above were obtained by Muzzio and coworkers |^p|] for a number of 
2D and 3D deterministic chaotic flows. Thus, the exact theory of the curvature statistics that we have been able to 
develop in the framework of the Kazantsev-Kraichnan model incorporates all of the essential features thus far observed 
numerically, as well as surmised in less direct theoretical ways. The close agreement between our theory and an array 
of numerical results obtained for more realistic flows provides an additional validation of our approach and, more 
importantly, suggests that the statistics of line- and surface-element advection possess a high degree of universality. 
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APPENDIX A: DERIVATION OF THE FOKKER-PLANCK EQUATION FOR THE JOINT PDF OF 

MAGNETIC FIELD AND LORENTZ TENSION 

Let us briefly describe the (standard) procedure we used to derive the Fokker-Planck equation (||). 
In the case of incompressible advecting flow, the magnetic field and the Lorentz tension satisfy 

O t B l +u k B l k =u l k B k , (Al) 
d t F l + u k F l k = u l k F k + u l km B k B m . (A2) 

We start by introducing the characteristic junction of the fields B(<,x) and F(i,x) at an arbitrary fixed point x, 

Z{t;n,\) = (Z(t,x;p,A)> = (cxpfz/i^t, x) + i\<F l (t, x)] ) (A3) 

Here and in what follows the angular brackets denote ensemble averages and overtildes designate unaveraged quantities. 
The function Z(t; /j,, A) is the Fourier transform of the joint PDF of the vector elements B l (t, x) and F l (t, x). Clearly, 
Z cannot have any spatial dependence due to the homogeneity of the problem. 

Upon taking the time derivative of the unaveraged function Z(t, x;/U, A) and making use of the evolution equa- 



tions (Al) and (A2), we find that Z satisfies 



i),Z + u k Z M = ( Mi A + u] k Z - *A 4 -^— u] km Z. (A4) 
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In order to establish an evolution equation for the (averaged) characteristic function Z(t;(j,, A), we must average 
the three mixed products of Z and the velocity field that appear in the above equation. The average that arises 
from the convective term vanishes due to the incompressibility of the velocity field and the homogeneity of the 
problem: (u Z t k) — ~( U \Z^ = 0. The remaining two averages are computed via the standard Gaussian splitting 
mechanism 

{<k z) - -<yz,) .-| tf /A>'(uK(f,l)A {awtf) = -5"& («£ + z - (A5) 

< A J) - <.%_> - jT ■* / dVMt xK( f <, x',)^ (^j) - - i <4U*^fc* CM) 

where we abbreviate k*£j = «*jy(y — 0), and n kmln = K km i n (y — 0). The above expressions have been obtained as 
follows. The functional derivative that appears under the integrals is the first-order averaged response function. It 
satisfies the causality constraint in that it vanishes for t' > t, whence follows the upper limit of the time integrations. 
Since the velocity field u l is 6 correlated in time, the time integration is removed and only the equal-time value of 



the response function has to be calculated. That is done by formally integrating Eq. (A4) from to t, taking the 
functional derivative 6/6vP(t',x!) of both sides, averaging, setting t = t' , and taking causality into account. The 
result is 



5Z(t,x) 



\5vP{t 

After integration by parts, the spatial integrations are removed due to the presence of 5 functions. Note that we make 
use of the fact that odd derivatives of the velocity correlation tensor K lJ (y) vanish at y = 0. 



Upon averaging both sides of Eq. (A4) and using the expressions (A5) and (A6) for the mixed averages, we obtain 



a closed evolution equation for the characteristic function Z, 

dtZ = ~ «« (fii^- + Xi^-] (m^- +X j ^ r )z-l K% lmn XiX j a * Z. (A8) 

2 M \ dfi k dX k J V dfn dXi) 2 - klmn dfi k dfj,idfj, m dfi n 

Inverse Fourier transforming this equation yields the desired Fokker-Planck equation (Q) for the joint one-point 
probability density function of the magnetic field B and the Lorentz tension F. 

The derivation of all other Fokker-Planck equations that appear in this paper follows the same general outline. 

APPENDIX B: COMPRESSIBILITY EFFECTS 

Let us relax the incompressibility condition and allow the advecting velocity field to possess an arbitrary degree 
of compressibility. Mathematically this means that we have to retain the terms involving divergences of u in the 
equations (Q) and (j^) and to allow the compressibility parameters a and b in the small-scale expansion (^) of the 
velocity correlator to vary in the intervals 

1 2 

<a<l, <6<2, (Bl) 

d+l - - ' d+3 ~ ~ ' v ; 

where the lower bounds correspond to the incompressible and the upper to the irrotational case. We will often use 
an alternative pair of compressibility parameters (3 = d[l + (d + l)a] and ( = d[2 + (d + 3)b] that have the advantage 
of being always nonnegative and vanishing in the case of incompressible velocity field. 

The exact treatment of the joint probability distribution of F and B is completely analogous to that presented 
in Sec. |l| for the incompressible case. The Fokker-Planck equation is now 

ap - -- K lj ( -5 k + —B k - 5 k —B r + —F k - 2S k —F 1 
0t ~ 2 K > kl { ° l + dB> 1 dBr + dF i* Zd * dF r t 

9 7 B l - 5'J-B S + J-F l - 2SiJ-F s I P 



dB'J ' Oil- dFi ■' ill- 
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The quantities (F 2 ) and (B 4 ) again satisfy equations (J?]) and (^J), respectively, with coefficients jf, Sf, and 74 
modified to include the dependence on the compressibility parameters (3 and £. The general expressions for these 
coefficients are listed in Table |. Note that the source term Sf remains positive for all allowed values of b. The 
steady-state solution of the form (^) continues to exist provided 7f — 74 < 0, which is satisfied for values of the 
compressibility parameter a such that 

d-2 

a<flc= 2(3^2) (B3) 

(in d = 2, a c = 0, in d = 3, a c = 1/14; this inequality can also be derived from a generalization of the simple 
argument in support of folding given in the Introduction: see Appendix ^|). Thus, for "nearly incompressible" flows, 
the folding picture persists in the strong sense that the parallel scale of the field remains approximately constant and 
comparable to the characteristic scale of the advecting flow. On the other hand, if the flow possesses a fair degree of 
compressibility, the parallel scales will start decreasing exponentially. 

Let us now retrace the path taken in Sec. |ll] and study the evolution of mean-square curvature and mirror force 
in the case of arbitrary degree of compressibility. Again, equations ( |li"| ) and ( |l2| ) preserve their form with modified 
coefficients jk, &k, 1m, 7m k, Sm (see Table |). None of these quantities changes its sign for any allowed values of 
the compressibility parameters. The essential structure of the solutions, therefore, does not change compared to the 
incompressible case, and the growing mean-square curvature (K 2 ) remains the one interesting quantity to watch. 

As we discovered from the statistics of the Lorentz tension, for a < a c the anticorrelation between the magnetic-field 
strength and the field-line curvature is preserved: while (F 2 ) / (B 4 ) remains constant, (F 2 /B 4 ) grows at the rate tr-. 
However, once the compressibility parameter a exceeds the critical value a c , the ratio (F 2 )/(B 4 ) starts growing as 
well, and the anticorrelation between B and K is weakened. Comparing the growth rate 7^- of the mean-square 
curvature with the growth rate jp — 74 of the ratio (F 2 )/(B 4 ), we find that jk > If — 74 provided 

3 

a < a * = 4d^7' (B4) 

While in 2D the second critical value a* = 3 lies outside of the interval of allowed values of a e [—1/3, 1], in 3D we 
have a* = 3/5 < 1, which is permitted. Thus, in three dimensions, for a > 3/5, the negative correlation between the 
field strength and the field-line curvature is replaced by a positive one, so the regions of maximal growth of the field 
and its curvature coincide! 



To prevent any misconception from arising with regard to the quantitative character of the conditions ( B3 ) and ( B4 ) , 
we ought to remark here that the particular critical values of the compressibility parameter when one or other statistical 
correlation breaks down are, of course, largely functions of what particular statistical averages are used to measure 
these correlations. Such sensitivity is due to the high degree of intermittency of the statistics of passively advected 
fields. 

Let us discuss the implications of the new facts that have emerged from this excursion beyond the confines of the 
incompressible advection theory. Clearly, the main feature of the compressible regime is that the velocity field is freed 
from having to preserve the volume and, along with stretching vortical motions that characterized the incompressible 
case, there now are motions that contract (or inflate) the volumes, with magnetic-field lines trapped inside. The struc- 
ture of the magnetic field now depends on the competition of stretching and contraction, whose relative importance 
is measured by the compressibility parameter a. We have seen in this section that stretching wins as long as a stays 
below a certain critical value a c . Once this value is exceeded, the parallel scale of the field cannot be prevented from 
decaying exponentially. While it may still be decaying slower than the perpendicular scale, thus giving rise to "small 
folds," both scales are now deep in the subviscous range and will eventually equalize when the resistive cut-off scale 
is reached. A tangled state will result. As a increases, the anticorrelation between the strength of the field and its 
curvature gradually weakens and, in 3D, is even reversed when a reaches a second critical value a*. This gives another 
indication of the increasingly tangled nature of the growing magnetic field. 



However, as is seen in Sec. [II, the growth of the magnetic field in sufficiently compressible flows only takes place in 
a small fraction of the total volume of the system, while elsewhere both the field strength and the field-line curvature 
remain relatively low. Thus, the tangled state is not set up everywhere throughout the system, but only in a small 
part of it where there is an appreciable growing magnetic field. This situation is, of course, due to volume contraction. 
The distribution of the density of the advecting medium is lognormal (highly intermittent) [ p2[ . While (p 2 ) and all 
higher density moments grow exponentially [Eq. ( p8|) of Appendix |c| , the growth of the density only occurs in a 
small fraction of the volume of the system. This is very natural and could not have been otherwise, for, as the total 
mass of the medium is conserved, (p) = const, the exponential growth of the higher moments of the density must be 
compensated for by the exponential contraction of the regions that are responsible for this growth. (In this context, 
one may also recall the results of Chertkov et al. US who found that for compressible-enough advecting velocity fields, 
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the Liapunov exponents for the Lagrangian fluid-particle separation become negative, so the fluid-particle trajectories 
tend to converge.) The density statistics are known to be intimately related to the statistics of the magnetic field 
(see Ref. |3^] and Appendix |(]) . Namely, there is a positive correlation between the density of the medium and the 
strength of the frozen-in magnetic field. This positive correlation can be deduced from the fact established in Ref. |[|] 
that even moments of are universal functions independent of the density statistics. The magnetic field will, 

therefore, tend to grow wherever the density does. 

Finally, let us note that, in the case of compressible MHD turbulence, the condition for the onset of nonlinearity 
is determined not just by the magnitude of the Lorentz tension B • VB, but rather by that of the total Lorentz 
force divided by the density of the fluid: (— X7B 2 /2 + B • VB)/p, which includes the magnetic pressure term. The 
latter will grow exponentially due both to the amplification of the magnetic energy and to the decrease of the (total) 
characteristic scale of the magnetic field. It is not hard to see that it will quickly outgrow the tension term and the 
nonlinear effect associated with magnetic pressure will dominate. 



APPENDIX C: AN ALTERNATIVE DERIVATION OF THE CRITICAL DEGREE OF 

COMPRESSIBILITY 



Let us demonstrate how the critical value of the compressibility parameter a derived in Appendix 
by constructing a generalization of the simple argument we gave in the Introduction [formula (^) 
continuity equation for the density p(t, x) of the medium, 



can be obtained 
Upon using the 



± p =-pV-u, (CI) 



we find that the magnetic field B(t, x) and the Lorentz tension F(i,x) satisfy 



12 = 2- Vu, (C2) 
at p p 



d F F BB BB 

- T = ^- Vu+ : VVu VV • u. C3) 

at p z p l p p p p 



We again suppose that the parallel variation of B is initially on scales much smaller than those of the velocity field u. 
Then the terms in Eq. ( |C3| ) that contain second-order derivatives of u are subdominant and can be neglected. We see 
that in such a case B/p and F / p 2 satisfy the same equation, which is the equation for the advection of a (contravariant) 
passive vector W, 

— W = W • Vu. (C4) 
at 

The statistics of passive vectors were treated (as a particular case of the statistics of general tensor fields) in Ref. jj2) . 
It was proved there that these statistics could be separated into two independent parts: one universal, the other nonuni- 
versal, the latter being expressible in terms of the statistics of the density. Specifically, the even moments of Wp 1 ^ 
are universal functions independent of the statistics of the density p. Therefore, the even moments of B/p 1 ^ 1 /^ and 
of F/p 2 ^ 1 ^ should also be independent of the statistics of the density. Making use of these results, one can write 

{B 2n ) = V) = const / d (M)(p 2n(d - 1)/d >, (C5) 

(F 2n ) = ((^) 2n p 4n ) = const / d (n,0<p 2 " (2rf - 1)/d ), (C6) 

where fd(n, t) = ((Wp 1 ^ d ) 2n ) is a universal function. Both fd(n, t) and the moments of the density field were calculated 
in Ref. 1331 , 



fd(n, t) = const exp 



1 n(2n + d)(l + a)n 2 t 



(C7) 



d 

(p 2n ) __ const exp [n(2n - 1)0 Kit] , (C8) 

where = d[l+a(d+l)} (vanishes in the incompressible case). With the aid of the above formulas, it is straightforward 
to calculate 
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=7||* 



7|| 



[2(3d-2)a-(d-2)]/s 2 . 



We see that 711 < for values of the compressibility parameter a such that 

d-2 



a < a r 



2(3d-2)' 



(C9) 



(CIO) 



We have thus recovered the inequality ( B2 ) . 

We would like to emphasize that the above derivation clearly demonstrates that the effects of compressibility on 
the field structure are due to the crucial part that the density of the medium plays in determining the statistics of 
the magnetic field in compressible flows. 
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